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main result is applicable to a general class of adiabatic responses. 



f~i ' Abstract. As a first step in the project of ratchet optimisation, the optimal driving 

waveform among a wide class of admissible functions for an overdamped, adiabatic 
C/3 . rocking ratchet is shown to be dichotomous. 'Optimum' is defined as that which 

I achieves the maximum (or minimum negative) average particle velocity. Implications 

'"^ ' for the design of ratchets, for example in nanotechnological transport, may follow. The 

O 

p • 1. Introduction 

Though its history stretches back to Feynman yj and Smoluchowski [2j, the ratchet 
^ I concept — a non-zero particle current arising due to symmetry breaking [3l H] of a 
potential or force — has experienced a resurgence of interest in the last 15 years. Initially, 
the interest was motivated by the possibility of the ratchet concept explaining the 
operation of biological molecular motors O O l^. Ratchets have now been realised 



O I in a wide range of physical systems [8j , with possible applications in nanotechnology [9] , 
I and are being explored in quantum [10] as well as classical regimes. 

Most authors choose a potential and a force that preserve or break the relevant 
! symmetries, and then investigate the dependence of the ratchet current on scalar 
^ I parameters, such as damping coefficient or temperature. A few have considered 
a parametrised potential or force, such as a sawtooth potential with variable up- 
slope and down-slope gradients [11], and then vary those parameters. Another such 
parametrisation is of the amplitudes and phase of a biharmonic force, as will be discussed 
in section m 

In this article I begin the ambitious project of optimising the 'shape', or functional 
form, of the symmetry-breaking force or potential over all admissible functions, to 
achieve maximum ratchet current. Previous discussions of ratchet optimisation have 
generally been in an abstract thermodynamic sense, pointing out that maximum 
efficiency of a Brownian ratchet is achieved in the limit of tight coupling [12i| . 
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It is unlikely that one shape will be optimal for all regimes: overdamped, dissipative 
and Hamiltonian; fast and adiabatic driving. One specific regime is here considered, this 
being an overdamped, adiabatic rocking ratchet, with equation of motion 

rjx{t) = -V'{x)+f{t)+at). 

This models an overdamped particle with position x{t) and damping coefficient t] moving 
in a spatially- varying potential V{x), subject to a time- varying force f{t) and Brownian 
noise C,{t). In the adiabatic limit the average velocity of the particle in the ratchet is [8] 



(i) 
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= rL , r.^L , . ... . \ TTTT^' (2) 



vIodxJ^^^dzexp{[V{z) -V{x) - {z - x)y\ /ksT} 
where the periods of the force and potential are Ty and L, respectively, and the 
temperature of the Brownian noise is ksT. 'Adiabatic' in this context means that 
the frequency of the driving is slow compared to any other characteristic frequencies of 
the system, such as the motion of the particle in its potential V{x). 

The potential is here set to the symmetrical V{x) = l^cos(27ra;/L), so that 
an average particle velocity arises from temporal symmetry breaking only, and the 
remaining parameters are scaled so that V = 1 and L = 2tt. In all calculations in 
this document rj = 1\ any other value, as can be seen from ([2]), will just scale the 
velocity (x). Plots of the resulting function v{y) are shown in Figure [TJ Notice that 
at large ?/, v{y) — > ?/; differences are only at small and that even these variations 
disappear as fc^T becomes large. 

The optimal force, as will be proved in section 2, is dichotomous, that is, it takes 
only two values, switching instantaneously between them. This result holds for the 
optimisation of any function v{y) for the measure ([1]) over the admissibility criteria 
presented in that section. Using some more specific properties of the response function 
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([2]), section 3 finds further properties of the optimal force for that response function. 
Section 4 presents numerical illustrations and general discussion, followed by conclusions 
and directions for future work in section 5. 

To my knowledge, only one author has previously considered the dependence of 
the ratchet current on the general shape of the potential or force. In a recent article, 
Chacon |13] admirably proposed a 'degree of symmetry breaking' (DSB) for each of 
the relevant symmetries, in an effort to quantify the connection between symmetry 
breaking and ratchet current. The only evidence he showed to support this measure, 
however, was an interpretation that for a biharmonic driving force it predicted the same 
optimal combination of amplitudes as the perturbation analysis discussed in section HI 
Furthermore, the measures diverge for a wide range of parameters, while any DSB, if it 
is to be correlated with ratchet current, should presumably remain finite. (The Cauchy 
principal value of the DSBs do exist, but, it can be shown, at least one is constant for 
all biharmonic forces, which is not much use.) DSBs will not be considered further here, 
except in the sense that the adiabatic response ([IH2]) provides an indirect measure of 
symmetry breaking, for the current is zero when all symmetries are unbroken. 

2. The optimal force is dichotomous 

The class of admissible functions f{t) over which the optimisation will take place is the 
set of (i) piecewise continuous functions with a finite number of discontinuities ('almost 
continuous'), which (ii) have period T, that is, f{t + T) = f{t), are (iii) bounded, that is, 
there exists an L for which |/(t)| < L for all t, and (iv) have zero mean, /J" f{t)dt = 0. 

Theorem 1. If force (s) exist that optimise the current 



then one such optimal force is dichotomous, that is, for all t, f{t) = M or N, for some 
M and N. The adiabatic response function v{y) is assumed to have u(0) = 0. 

An 'optimum' force is defined as a force f* with {x){f*) > {x){f) for maximisation 
or {x){f*) < {x){f) for minimisation, where f is any admissible force. 

It is trivial to see that the adiabatic response function ([2]) satisfies f (0) = 0. 

Proof. First consider a constant force. This, by the constraint of zero mean, must have 
value zero, which gives zero velocity, and is therefore not optimal. Therefore any optimal 
force must take two (dichotomous) or more values. 

Choose a small discretisation size dt to discretise the force f{t) into a series of 
segments f{[n + ^]dt) on ndt < t < {n + l)dt, n = 0, . . . ,T/dt — 1. Let /, refer to the 
discretisation of f{t) at tj. Consider an admissible function which is not constant and 
not dichotomous, that is, it takes three or more values. For a small enough dt there 
must exist three tj, i = 1, 2, 3 with /j ^ fj for i ^ j, that is, the force takes on at least 
three different values. The contribution of the three intervals to the velocity is 



(i)(/) ^ / vifit))dt 
Jo 




(3) 
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Figure 2. Schematic of tlie metliod to convert a trichotomous region of tlie discretised 
force into a dichotomous region. Having chosen the three values /i,2.3 in (a), the values 
of the force at those times are replaced with just /i and /2 in (b). In total, the force 
takes values /i_2 for total times dt'i 2, respectively. 



Consider the effect on the velocity of replacing this three-valued section with just 
two values. Suppose /s is to be eliminated. This procedure is illustrated in Figure Hi 
Let the combined lengths of the intervals at which f{t) will be set to /i and /2 be dt[ 
and dt2, respectively, where dt'i + dt2 = 3dt. For Jq f(t)dt to remain at zero, 

(/i + /2 + fz)dt = f\dt\ + /2(i4. (4) 

At least one fi must be greater than the average value (/i + /2 + /3)/3 and at least 
one less than the average value. If, without loss of generality, these are labelled as /i 
and /2, respectively, then a solution dt'^, dt'2 to (jl]) can always be found. Using (jl]), the 
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Figure 3. Cases to be considered for tlie configuration of tlic points {fi,v(fi)),i = 
1,2,3. (a) Tlie three points are collinear. (b) The point with intermediate y-value 
(note y is the horizontal co-ordinate here) is below the line joining the other two 
points, (c) The point with intermediate y-value is above the line joining the other two 
points, and has y- value less than the average value (/i + /2 + /3)/3. (d) The point with 
intermediate y-value is above the line joining the other two points, and has y-value 
greater than or equal to the average value /av = (/i + /2 + fs)/^- The points have 
been numbered so that /i > /av, /2 < /av, and removing in the procedure outlined 
in Figure [2] will increase the ratchet velocity. 



change in velocity due to the change from a trichotomous to a dichotomous region can 
be shown to be 

A(x) = (rft/T)(/3 - /2)(m2i - m23), (5) 

where 

- /, - /. 

is the gradient of the hne joining {fi,v{fi)) and {fj,v{fj)). 

Each of the four possible arrangements of the points v{fi)) illustrated in Figure 
[3] will now be considered in turn. 

Case illustrated by Figure\3(a). In the collinear case there will be no change in 
velocity, A(i;) = 0. 

Case illustrated by Figure \B(b). Without loss of generality the points may be 
labelled as shown. The middle point will therefore be 'removed'. As can be seen on the 
figure, 77121 > "^23, and since /a > /2 then from ([5]), A(i;) > 0. 
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Case illustrated by Figure \^c). Without loss of generality the points may be 
labelled as shown, that is, the most negative / as '3'; this is the point that will be 
'removed'. From the figure, m2i < ^23, and since /s < /2 then from ([5]), A(x) > 0. 

Case illustrated by Figure \3(d). Without loss of generality the most positive / 
may be labelled as '3'. From the figure, m2i > m23, and since /s > /2 it follows that 
A{x) > 0. 

Cases (c) and (d) are treated differently so that in each there exists a solution to 
(jlj), as discussed in the text below that equation. 

Thus replacing the trichotomous (three-valued) region by a dichotomous one will 
increase the ratchet velocity or at worst keep it the same. In doing so the A^-valued 
force has been replaced by an (A^ — l)-valued force. Repetition of the above procedure, 
at each step increasing the ratchet velocity, will end with a dichotomous force. (Some of 
these iterations will require operation on segments with unequal width, but the above 
mathematics can be easily modified to account for this.) As the width of the original 
discretisation dt approaches zero, the above method approaches an optimisation of the 
original waveform. 

For any admissible force, then, a dichotomous one with higher velocity can be found. 
Therefore an optimal force for maximum velocity, if it exists, must be dichotomous. By 
an analogous argument to the above, another dichotomous force must provide the most 
negative (minimum) velocity. □ 

The procedure in this proof will not necessarily converge to the optimum 
dichotomous waveform, but does always find one that gives a greater velocity than 
the arbitrary admissible waveform with which it began. 



3. Dichotomous optimal force for the ratchet 

This section finds further properties of the dichotomous, optimal force, assuming that 

d V (11 ) 

sign?/ > for ?/ 7^ 0, (6) 

dy y 

and that v{y) is continuous and smooth. The response ([2]) satisfies these conditions. 

The defining characteristics of a dichotomous force are the force's two values and 
the time spent at each value. The ratchet velocity (x)(/) is clearly insensitive to time 
shifts. Likewise it is insensitive to 'time-mixing' of the force, where it alternates many 
times between its values, as long as the total time spent at each value is the same, 
though a force with just two transitions per period is more physically realistic. 
Define the dichotomous force 

_fM>0, 0<t<Td 
•^^ '* ~ j < 0, Td<t<T ' 

Since the force must have zero mean, 

Md + N{1 -d) = 0. (7) 
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The velocity can therefore be written as 

(*)dichot=^W^ + ^(-^ 

Setting = ^j-, treating M as constant and using constraint 
function of d gives 

> M-N- • 

where the prime here denotes a derivative. This provides a condition for the locally 
extremising A^, A^*. It holds for all response functions for which the optimal force is 
dichotomous. It can be shown that oc v"{N), so the type of the extremum is 

determined by the sign of v"{N*). 

The next step is to attempt optimisation with respect to M. Return to ([HD. When 
forming one must consider dd/dM; however, these terms cancel to give 



(8) 

to specify iV as a 
(9) 



d{x) 
'dM 



d [v\M) - v\N*)] 



(10) 



This quantity, as will now be shown, is greater than zero for all M > 0. Notice 
first that ([9]) can be interpreted graphically: The tangent of v{y) at A^* is collinear with 
the secant joining {N* ,v{N*)) and {M,v{M)), as sketched in Figure HI An equation 
describing this line, to be denoted by T, is 

VTiyr) = v'{N*){yT - A^*) + v{N*). 

It intercepts the vertical (v) axis at 

vt{0) = -v'{N*)N* + v{N*) > 0, 

with the inequality following from ([6]) and A^* < 0. Since the line also contains 
{M,v{M)) then another equation describing it is 

VT{yT)=v'{N*){yT-M)+v{M). 
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According to this equation the f-intercept is 

vt{0) = -v'{N*)M + v{M). 
Suppose v'{M) < v'{N*). Then 

vt{0) < -v\M)M + v{M) < 

with the last inequahty following from ([6]) and M > 0. But this contradicts the earlier 
result. Therefore it follows that in fact v'{M) > v'{N*) and by that > for all 
M > 0. Therefore the optimum M, for maximum velocity, is its maximum permitted 
value, L. 

In order to find the waveform that minimises the velocity the above argument may 
be repeated but first optimising with respect to M. It then follows that the optimum 
is the most negative possible, —L. If v{y) = —v{—y), as does ([2]), then this force is 
the negative of the one just found to maximise the velocity. 

By inspection of Figure [H, it appears that for each v{y) given by ([2]) there exists 
exactly one A^* that satisfies ([9]) for fixed M, as sketched in Figure HI Further, it appears 
that v"{N*) is of the correct sign for local maximisation (in Figure Hj) or minimisation. I 
state without proof that this is indeed the case. Since there is only one local extremum 
in ratchet velocity (as one varies A^) and it is of the correct type, and v{y) is smooth 
and continuous, it follows that the global extremum is also here. 

By Theorem 1, then, there exists an optimum force for maximisation (and another 
for minimisation), it is dichotomous, and its characteristics are unique. 

4. Numerical work and Discussion 

Figure [5] shows the optimal A^* and d at various force limits L and temperatures ksT. 
Notice that, as intuitively expected from Figure HI for most L at /c^T = 0.1 the optimal 
A^* is on or close to the knee in Figure [Tj while at higher temperature A^* changes 
more with L. As the upper limit increases, the duty cycle shows that the dichotomous, 
optimal force spends a smaller fraction of time at this upper value. 
A common ratchet driving waveform is the biharmonic drive, 

•^bihar(^) ~ ^ cos tut + B cos{2Ljt + (p). 

According to perturbation analysis this will generate, in an overdamped spatially- 
symmetric ratchet, at the limit of small force amplitude, an average ratchet velocity 
that scales with A^Bcoscj) [l3l[T3l[3]. The optimal biharmonic phase is therefore = 
(or = TT, for transport in the negative direction) and, if the maximum value of the 
force is limited to L, the optimal amplitudes are B = L/3 and A = 2L/3. 

The velocities predicted by ([l]l2l) are shown in Figure [6] for both the optimal 
dichotomous force obtained previously and this perturbatively optimal biharmonic force. 
At small L the biharmonic responses are almost as large as the dichotomous responses, 
as is appropriate for a perturbation analysis. For larger forces, however, the biharmonic 
responses are significantly less than the corresponding dichotomous responses, and even 
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Figure 5. Optimal value for the negative part of the dichotomous force (when 
maximising the ratchet velocity), iV* (solid lines), and optimal value for the 
dichotomous force's duty cycle, d (dashed lines), as functions of the limit L, for two 
temperatures ksT. 
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Figure 6. Velocities, normalised by the force's period T/, predicted by (IT][2]) for the 
dichotomous optimal force and the corresponding perturbatively optimal biharmonic 
force with same maximum value L, for temperatures fcsT = 0.1 (solid lines) and 1 
(dashed hnes). 



decay for very large forces. While the rates of increase of the dichotomous responses do 
slow, they will show no such turning point, as shown in the paragraph following (JTOl) . 

For any experimentalist or engineer designing an overdamped ratchet close to the 
adiabatic limit — or even as a first guess away from the adiabatic limit — and who can 
freely choose the time dependence of the force up to a limiting value ±L, a dichotomous 
force, it has been shown, will give the largest possible average ratchet velocity. The 
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optimal characteristics of the force can be solved with the equations or approximated 
with the graphs given here. 

That the optimal force is dichotomous is to be intuitively expected when one notes 
that the response ([1]), as should be expected for an adiabatic response, depends only on 
the force /, not any of its derivatives or ant i- derivatives. There is no 'penalty', then, 
for the discontinuities in the dichotomous force. 

It is acknowledged that practical restrictions of specific experimental arrangements 
may in some cases more strongly influence the form of the force, for example biharmonic 
forces from the mixing of laser beams [15], or near-dichotomous forces from the 
application of pressure across a membrane [9] . 



5. Conclusions and future work 



An optimum force for an overdamped, adiabatic ratchet, amongst the admissible 
forces — briefly, forces with zero mean and maximum absolute value not greater than 
L — if it exists, is dichotomous. The optimum characteristics of the dichotomous force, 
for maximum velocity, are (any time-shifted or time-mixed version of this force will give 
the same velocity) 

' L, 0<t<d 
^ N*, d<t<T 

where N* < is a solution of 

^ <L) - viN') 

^ ^ L- N* 

and d is then given hj d = N* / {N* — L). From plots of the response ([2]), it appears 
that there a (unique) dichotomous optimal force does exist, though this was not formally 
proven. 

The main result of this article. Theorem 1, may apply to a wide range of problems. 
It was shown that for any response of the form Jq v{f{t))dt, where f{t) has a finite 
number of discontinuities, has zero mean and has size constrained by its maximum 
absolute value, the response is extremised by a dichotomous force. 

The result may be useful for engineers designing overdamped Brownian motors 
in the adiabatic limit, or even as a first guess away from the adiabatic limit, wishing 
to achieve the maximum possible average velocity from the ratchet. For large forces 
a dichotomous force can generate a much larger velocity than that generated by the 
biharmonic force optimised by perturbation analysis. 

There remain many ratchet configurations in which this optimisation task could be 
addressed. For the case of an overdamped, adiabatic ratchet with time-symmetric force, 
early results indicate that the optimum potential may be close to a sawtooth waveform, 
which would be a pleasingly complementary result to the above. One could consider an 
overdamped ratchet with fast driving, from which some conclusions about the optimal 
waveform for all driving frequencies might be drawn. There are also the underdamped. 
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both adiabatic and fast-driven, cases. In these analyses, more sophisticated methods 
such as the calculus of variations or Pontryagin's Maximum Principle [16] may be useful. 
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